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secondary antiproton and diffuse gamma-ray yields. 
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Recent data from CREAM seem to confirm early suggestions that primary cosmic ray (CR) 
spectra at few TeV/nucleon are harder than in the 10-100 GeV range. Also, helium and heavier 
nuclei spectra appear systematically harder than the proton fluxes at corresponding energies. We 
note here that if the measurements reflect intrinsic features in the interstellar fluxes (as opposed to 
local effects) appreciable modifications are expected in the sub-TeV range for the secondary yields, 
such as antiprotons and diffuse gamma-rays. Presently, the ignorance on the origin of the features 
represents a systematic error in the extraction of astrophysical parameters as well as for background 
estimates for indirect dark matter searches. We find that the spectral modifications are appreciable 
above 100 GeV, and can be responsible for ~30% effects for antiprotons at energies close to 1 TeV 
or for gamma's at energies close to 300 GeV, compared to currently considered predictions based on 
simple extrapolation of input fluxes from low energy data. Alternatively, if the feature originates 
from local sources, uncorrelated spectral changes might show up in antiproton and high-energy 
gamma-rays, with the latter ones likely dependent from the line-of-sight. 



PACS numbers: 98.70.Sa 

I. INTRODUCTION 

A more accurate determination of primary cosmic ray 
spectra at the top of the atmosphere has obvious im- 
plications for the understanding of the acceleration and 
propagation of Galactic cosmic rays. It is also crucial 
for other fields of investigations in astroparticle physics, 
two notable examples being atmospheric neutrino stud- 
ies (e.g. [![) and the calculation of the backgrounds for 
indirect dark matter (DM) searches (see for example 0). 

In the specific case of indirect DM searches, an impor- 
tant implicit assumption is that fluxes measured at the 
top of the atmosphere, at sufficiently high energies to 
avoid solar modulation effects, are representative of in- 
terstellar medium (ISM) spectra. Or, more correctly, one 
often assumes universality for the injection term and the 
propagation properties, fitting the free parameters (like 
injection and diffusion power-law index) in such a way to 
reproduce the observed spectra. It is those "universal" 
interstellar spectra which in turn enter as source term of 
secondary yields (like antiprotons or diffuse gamma rays) 
due to inelastic collisions in the ISM. This assumption is 
usually supported by the apparent featureless nature of 
the observed cosmic ray fluxes (suggesting, at least in 
average, some universal mechanism for production and 
propagation) as well as by the check a posteriori that 
the diffuse gamma-ray radiation of hadronic origin has 
a spectrum consistent with the hypothesis, within the 
errors. However, especially at energies larger than the 
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TeV scale, inferring accurate spectra is challenging due 
the scarce statistics and experimental difficulties, mak- 
ing the above arguments at best based on shaky observa- 
tional evidence. Also, for assessing uncertainties in sec- 
ondary yields, the usual practice is to fit primary data 
to some power-law parameterization and extrapolate to 
high energies. While this is a reasonable prescription 
for most applications given the present level of under- 
standing, these simplified approaches and assumptions 
might hide a systematic error when searching for signa- 
tures showing peculiar energy features. For example, for 
antiprotons this is the case involving contributions from 
DM annihilation Q or production at the sources Q. 

Obviously, the standard prescriptions do not usually 
account for the possibility that a systematic departure 
(rather than statistical scattering) is present in the spec- 
tral shape of the fitting formula, which is mostly cali- 
brated on low energy data, neither of the possibility that 
observed spectra might not be fair representative of the 
interstellar ones. Recent data from the CREAM balloon- 
borne experiment [f| seem to confirm earlier suggestions 
(see e.g. [6|]) that cosmic ray spectra at few TeV/nucleon 
are harder than in the 10-100 GeV range, and that he- 
lium (He) and heavier nuclei fluxes are harder than the 
proton (p) flux at corresponding energies. Preliminary 
data from PAMELA also suggest a hardening in p and 
He spectra at a rigidity of about 250 GV, with a He spec- 
trum having an index ^0.1 lower than the proton one 
over all energies above a few GeV [3] ■ All this motivated 
us to have a more careful look at the errors potentially 
committed when estimating secondary yields in the in- 
terstellar medium. 

In this article, we refrain from discussing possible as- 
trophysical interpretations of the above mentioned fea- 
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tures, although some have been proposed, see d, We 
note however that if the measurements reflect intrinsic 
properties of the interstellar spectra, appreciable mod- 
ifications (i.e. above ~ 10%) of very specific spectral 
shape are expected for the secondary yields in the 0.1 to 
1 TeV range, which is directly accessible (with growing 
precision) to present and forthcoming experiments like 
PAMELA, FERMI, and AMS-02. On the other hand, 
if the hardenings reflect local phenomena/sources, sec- 
ondary yields which probe a large volume of the ISM, 
like antiprotons, might not show relevant departures from 
naive expectations, while the diffuse gamma-rays along 
different lines of sight might reveal different hardenings 
reflecting the primary spectra present in different regions 
of the ISM. Clearly, this provides an important test for 
theories about the origin of the breaks. To the best of 
our knowledge, present data in high energy astrophysics 
are either unrelated to the hardenings discussed here or 
still of too limited precision to provide a crucial test, but 
the situation is likely to change in the near future. 

This article is structured as follows: in Sec. HJwe dis- 
cuss the input fluxes and parameterization used to pro- 
vide a first estimate of the effect. In Sec. Mil we present 
the results for antiprotons and 7-rays, finally in Sec. IIVI 
we discuss some implications of our findings, and con- 
clude. 



II. INPUT FLUXES 

In the present exploratory study, we refrain from the 
ambitious goal of analyzing the whole body of cosmic 
ray flux data in the 10 — 10 4 GeV/n range. Rather we 
limit ourselves to provide a first assessment of the sys- 
tematic effect potentially introduced by deviations from 
the power law behaviour at high energy, in general with 
different spectral indexes for different species. To this 
purpose, we explore the effects of combining the fits of 
"low-energy" (namely in the range about 10-100 GeV/n) 
proton {i — 1) and helium (i — 2) flux data, <pf, taken 
from AMS-01 0] (in turn, to large extent consistent with 
what reported by other experiments), with the "high- 
energy" (above about 1 TeV/n) fluxes 4>f inferred by 
CREAM @. We adopt broken power-laws to connect 
the two sets, using the following flux parameterizations 
(differential fluxes with respect to kinetic energy per nu- 
cleon T): 

MT) = ^(T)e(B 1 -T) + ^(T)Q(T-B 1 ), (1) 
cp 2 {T) = ^(T)e(B 2 -T) + ^(T)Q(T-B 2 ), (2) 

The fluxes "L" are the best fit values taken from AMS- 
01 @, rewritten in terms of kinetic energy T per nu- 
cleon (in GeV/n) instead of rigidity and asymptotically 
decreasing as ~ y-^ 'S f or ^ an( ^ y-2.74 j? or jj e j n un its 

of (GeV/nm 2 ssr) -1 , they write 

<t>i(T) = 1-71 x 10 4 ^(T + m p y-ml) , (3) 



^ (4r + W 2 Hc)2 - m S£j (4) 

The high energy fluxes "H" are taken from CREAM with 
the following criteria: i) power-laws in T are assumed, 
with the spectral indexes fixed to the best-fit values re- 
ported in |5j, i.e. 2.66 for p and 2.58 for He; ii) the proton 
spectrum normalization is taken from the first CREAM 
point in Fig. 3 of [5j; iii) the Helium spectrum normal- 
ization follows from imposing that at T — 9 TeV/nucleon 
the proton to helium flux ratio is equal to 8.9 [5]. As a 
result, in units of (GeV/nm 2 ssr) -1 , 

4>i{T) = 7.42 x 10 3 T~ 2 ' 66 , (5) 
4>%{T) = 4.03 x 10 2 T~ 2 ' 58 . (6) 

The crossover energies B\ , B 2 for the broken power-laws 
are simply obtained by continuity, and are approximately 
B p = 1000 GeV, B Hc = 30GeV/n for the parameters 
above 1 . A comparison with the predictions following 
from the extrapolation of the AMS-01 fits (i.e. the tfif of 
Ecis. (|3I4| )) to arbitrarily high energy will be presented to 
provide an estimate of the impact of high-energy spectral 
uncertainty on the secondary yield flux. 



III. RESULTS 

Discrepant hardenings of primary cosmic ray fluxes, 
possibly of non-universal nature, would obviously affect 
all the yields of e + , p and 7 secondaries produced by 
collisions in the interstellar medium (ISM). Here we do 
not discuss charged leptons simply because the primary 
flux effects do not provide the major uncertainty in the 
flux shape (even fixing the av erage p ropagation parame- 
ters): very likely recent data |KH12| indicate that addi- 
tional sources of "primary" positrons exist for which the 
above mentioned effects are expected to be sub-leading 
(see e.g. [l3[ ). Additionally, energy losses make the range 
shorter and the computation of the actual flux at the 
Earth non-trivial, so it would be more difficult to dis- 
entangle the effects due to the break in primary spectra 
from a complicated interplay of effects involving the dis- 
creteness of local sources, inhomogeneities in the radia- 
tion field, etc. as illustrated for instance in |l4j. The 
effect of universal primary CR hardening should be ap- 
preciable in the predicted shape of the antiproton or dif- 
fuse gamma-ray signal. Here we report a careful com- 
putation of the effect on the antiproton spectrum, where 
the impact is expected to be the largest in view of future 
high-statistics results from AMS-02, and an estimate of 



Assuming a relative uncertainty in the flux normalization of the 
two experiments of < 20% — certainly consistent with published 
values — would suffice to bring these crossover values in consis- 
tency with the rigidity ~ 250 GV hinted to by PAMELA, ref. 



3 




FIG. 1: Ratio of antiproton fluxes from hard sources (Eqs.JT] 
O) to the same flux obtained with p and He extrapolated 
from AMS data to all energies (see text for details) . 



the effect on the hadronic gamma-ray diffuse background, 
of some interest for the interpretation of FERMI data. 



A. Effects on secondary antiprotons. 

The computation of the secondary p flux has been per- 
formed as described in Refs. 0, H|, to which we refer 
for all the details. The only component which we will 
modify in the present calculation is the input p and He 
spectra. We briefly remind that secondary p are yielded 
by the spallation of cosmic ray proton and helium nu- 
clei over the H and He nuclei in the ISM, the contribu- 
tion of heavier nuclei being negligible. The framework 
used to calculate the antiproton flux is a two-zone dif- 
fusion model with convection and reacceleration, as well 
as spallations on the ISM, electromagnetic energy losses 
and the so-called tertiary component, corresponding to 
non-annihilating inelastic scatterings on the ISM. The 
relevant transport parameters are constrained from the 
boron-to-carbon (B/C) analysis [HI and correspond to: 
i) the half thickness of the diffusive halo of the Galaxy 
L; ii) the normalization of the diffusion coefficient Kq 
and its slope S (K(E) = K (3R S ); hi) the velocity of the 
constant wind directed perpendicular to the galactic disk 
V c = ±V c e* z ; and iv) the reacceleration intensity param- 
eterized by the the Alfvenic speed V a . The above pa- 
rameters show significant degeneracies when confronted 
to B/C data [15|]. Nevertheless, the impact on the sec- 
ondary p flux is marginal 0. The fluxes presented be- 
low have been obtained for the B/C best fit propagation 
parameters, i.e. L — 4 kpc, Kq — 0.0112 kpc 2 Myr _1 , 
6 = 0.7, V c = 12. km s" 1 and V a = 52.9 km s" 1 [ijf . 

We arc interested in the effect of primary p and He 
hardening at high energies on the p flux and therefore 
concentrate on the relative shape effect through antipro- 
ton flux ratios. Our results are reported in Fig. [TJ where 



we plot the ratio of antiproton fluxes obtained with two 
different primary spectra. The flux at the numerator has 
been obtained with the spectra in Eqs.([TJ[2]), while in the 
denominator we employ the fit to AMS data arbitrarily 
extrapolated to the highest energies. The modification 
of the antiproton flux clearly reflects in its shape. The 
effect of the hardening of primary spectra at hundreds 
of GeV/n starts to be visible on the antiproton flux at 
around 100 GeV. It is near 15% at 200 GeV and reaches 
30% at 1 TeV. Given the weak dependence of the sec- 
ondary antiproton flux on the B/C selected transport 
parameters, our results can be considered nearly inde- 
pendent of the propagation model. If the hardening of 
primary nuclei will be confirmed at high energies, a spec- 
tral distortion of the secondary antiproton flux has to 
be expected. This effect could be potentially observable 
by a future high precision space-based mission, such as 
AMS-02. 



B. Effects on hadronic diffuse gamma-rays. 

The cosmic gamma ray flux observed in our Galaxy 
is expected to be mainly due to the inelastic scattering 
of incoming CRs on the nuclei of the ISM. The involved 
hadronic reactions produce gamma rays mostly via ir 
decays. In addition to this hadronic component, other 
contributions are expected — at different levels depending 
on the specific model — to Inverse Compton and brems- 
strahlung radiation. The basic models for the production 
of gamma rays from ir° decays, considered for example 
by the Fermi-LAT Collaboration [l6j], do not introduce 
high-energy spectral breaks in the proton spectrum 0i, 
and account for nuclear effects (both in CR spectra and 
in target composition) in the ir° yield simply by reseat- 
ing the pp production via a constant "nuclear enhance- 
ment factor" , taken from the value at the reference en- 
ergy =10 GeV/n reported in [l7|. This enhancement 
encodes the relative yield of gamma-rays from nucleus-p 
and nucleus-Helium collisions compared with that from 



p-p collisions via appropriate factors , rrii a , basically 
constant at T > lOGeV/n (the effects discussed in this 
paper are only relevant at high energy, so it's enough to 
focus on quantities at T > 10 GeV/n). This enhancement 
is defined as 



mu 



1 



(7) 



where the index i runs over all CR species (including 
protons, i = 1), r ~ 0.096 is the He/H fraction in the 
ISM and 0i being the CR spectrum of the species i. If 
all the nuclei have roughly identical T dependence of their 
spectra, as suggested in [5| , one can write 



e M (T) = 1 



mi2 r 



1 



r 



. m 2 2r\ 2 , <t>N ,„s 

TO 2 i + )-r + k N——, (8) 

1 - r ) (pi 0i 



where 4>n(T) is any nuclear-like CR flux, and is a 
normalization factor. In Fig. [5] we show 6m(T) for three 
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FIG. 2: Enhancement factor, see Eqs. (|7I8[1 . f° r the three 
representative cases described in Sec. MI Bl 



cases: i) the constant value em = 1.84, adopted for ex- 
ample in [l6| (long-dashed, black); ii) the fluxes of p and 
He are set to the broken power-law functions described 
above, while the last term fc/v</>iv/0i is taken constant 
in energy and fixed so that €m(T*) = 1.84 (short-dashed, 
blue), iii) As in ii) for p and He, but assuming for nu- 
clei heavier than He a constant contribution to e below 
200GeV/n (so that e M (T*) = 1.84), then rising as T 01 , 
as suggested by CREAM data (solid, red). 

In Fig. [31 we show the result of computing the diffuse 
gamma ray spectrum (via the kernel provided in 18] 2 ) 
using the AMS-01 spectral fits (f>f, extrapolated to arbi- 
trarily high energy (long-dashed, black curve). The flux 
has been multiplied by E^' 78 to underline the departure 
from identical power-law behaviour between photons and 
parent CR due to production cross section/multiplicity 
effects. Instead, if one keeps £m = 1-84, but introduces 
the broken power-law spectrum for the protons only as 
from Eq. (fTJ, around 300 GeV one would obtain ~ 10% 
higher gamma fluxes, as shown by the long-dashed, pur- 
ple curve in Fig. [31 This case is introduced in order to 
gauge visually the effect of the break of 2.78 — 2.66 ~ 0.12 
in the spectral index, between AMS and CREAM de- 
termination of proton spectra (to be compared with the 
~ 0.01 and 0.02 fit errors, respectively, reported by the 
experiments). The solid, red curve shows the effect of 
"discrepant hardenings" of the spectra, namely the T 
dependence of 6m- This constitutes the major distortion 
and is mostly due to He (as shown by the short-dashed, 



FIG. 3: 7-ray spectrum: the standard departure from equality 
of power- law with parent flux (dot-dashed, black), of adding 
the hardening in the p spectrum at TeV scale as suggested 
by CREAM (long-dashed, purple), of assuming the CREAM 
hardening for both p and He (short-dashed, blue), and of in- 
cluding the small effect of other nuclei as well (solid, red). 



blue curve); overall, the spectrum around 300 GeV is 30% 
higher with respect to naive expectations. 

The effect discussed here is already "an estimate of er- 
ror": assessing the error on this quantity goes beyond 
our purpose. However, we can safely conclude that our 
results are not significantly affected by the statistical er- 
rors with which the normalizations or spectral indices are 
known. We have checked this explicitly as follows: while 
keeping the Helium fluxes at low and high energies at the 
values described above, we have varied the normalization 
of the p fluxes at low/high energy in such a way that the 
p to He ration varies within ±1 (i.e. twice the statistical 
error) of 18.8 at 100 GeV/nucleon H, and within ±0.6 
of 8.9 at 9 TeV/nucleon [5]. The resulting variations in 
the shape of secondary radiation are negligible (i.e. at 
most at a few percent level). This is due to the fact that 
the (relatively small) effect of p flux renormalization and 
the change in €m anti-correlate, and tend to cancel each 
other. On the other hand, the exact value of the spectral 
hardening is more important: Fig. [3] shows that more 
than 1 /3 of the hardening is due to the assumed "best 
fit" spectral index difference of ±0.12 between low and 
high energy. This should be compared with the statistical 
errors of about ~ ±0.01 and ±0.02 quoted by the AMS 
and CREAM collaborations, respectively. 



2 Note that we are only interested in the effects that different high- 
energy CR spectra have on the gamma-ray spectrum at E-y ^> 
1 GeV, so the simplified formalism presented in [l8| and valid in 
the high-energy regime is sufficient for our purposes. 



IV. DISCUSSION AND CONCLUSIONS 

In this article we have argued that departure at high 
energy from a simple and universal power-law for all cos- 
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mic ray spectra, as suggested by recent data, should cause 
a spectral distortion in the spectra of secondary cosmic 
ray yields (like diffuse photons and antiprotons) com- 
pared to the predictions obtained extrapolating the best 
fits to low-energy data sets. We have illustrated this ef- 
fect using the best fit results of AMS-01 data at low ener- 
gies and the CREAM data at high-energy, finding effects 
exceeding 10% above ^100 GeV, and reaching about 30% 
for photons around 300 GeV and for p close to TeV en- 
ergy; this figure is somewhat sensitive to the systematic 
error on the spectral index at high energy as well as other 
eventual systematics which do not cancel out in ratios of 
species (like p/He). If the hardening in nuclei data would 
be due to local effect and not representative of the ISM 
average, the effect on the antiproton yield may be small 
while the hadronic diffuse gamma-ray spectrum could be 
modified differently according to the line-of-sight. No- 
tice that the effect of a possible harder nuclear spectra 
on atmospheric neutrinos was already estimated in 

One might wonder how relevant is a high-energy effect 
of a few tens of percent in a field where data are usu- 
ally plagued by larger errors. We think that, at present, 
this level of accuracy is becoming crucial for at least a 
couple of reasons: Firs t, sp ace experiments like FERMI 
or the future AMS-02 hfl are introducing us to a new 
era of large exposures, which can reveal more subtle fea- 
tures than previous cosmic ray or gamma-ray experi- 
ments. Fermi data errors at E 1 ~ 100 GeV are already 
~ ±20% [13] , and forecasts that have been presented sug- 
gest that AMS-02 (if performing close to specifications) 
will be certainly sensitive to effects of this magnitude, see 
for example [2l|. Second, both diffuse gamma-rays [22[ 
and the combination of hadronic data [23] are consistent, 
at least at leading order, with a "standard" scenario for 
the production and propagation of cosmic rays in the 
Galaxy. It is very likely that any departure from base- 
line models, if detectable, is going to be present at such a 
sub-leading level. Modelling thus the astrophysical back- 
ground for indirect DM searches as a simple power law, 
as often done in the literature, might lead to wrong con- 
clusions about the evidence of a signal, or to a bias in 
the inferred values of the parameters describing the new 



phenomena, should they be detected. 

Even in a conservative scenario, the detection of such 
spectral signatures in secondary channels would provide 
a way to check the interstellar nature of the spectral fea- 
tures in the cosmic ray flux at the Earth suggested by the 
present experiments. We believe that secondaries provide 
an important handle for an empirical cross-check. One 
should also consider the partial degeneracy of such effects 
with the extraction of propagation parameters, in order 
to fully exploit the statistical power of forthcoming data 
sets. Knowing better the primary flux shapes would al- 
low one to set strategies minimizing these effects. Last 
but not least, a multi-messenger approach would allow 
one to disentangle these features from alternative sources 
of spectral distortions: features similar to the ones dis- 
cussed in this article arise e.g. in models where high en- 
ergy p are produced in sources [|| , but in that case also 
associated signatures in secondary/primary "metals" [24| 
(and possibly in high energy neutrinos |25j) are expected, 
which are absent for the process described here. 

While we are entering a much higher precision era 
in cosmic ray studies, it is important to keep in mind 
a couple of points: i) that multi-messenger and multi- 
channel analyses are mandatory, if one is to gain some 
deeper knowledge of cosmic ray astrophysics, ii) That 
any hope for the detection of new physics (not to speak 
of extracting new physics parameters) requires a more 
robust understanding of the possible range of astrophys- 
ical yields. In that respect, a natural development of 
this initial investigation would be to (re) assess how the 
errors on primary flux knowledge map into the predic- 
tions for secondaries (including their normalization), as 
much as possible in a parameterization-independent way. 
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